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Abstract 

Electromagnetic  material  characterization  is  the  process  of  determining  the  com¬ 
plex  permittivity  and  permeability  of  a  test  sample.  The  primary  goal  of  this  thesis 
is  to  develop  a  new  two  transmission  material  measurement  method  to  decrease  the 
error  associated  with  using  a  reflection  measurement.  The  transmission  method  uses 
a  sample  transmission  measurement  and  a  acrylic  backed  sample  transmission  mea¬ 
surement.  This  technique  is  first  demonstrated  in  a  rectangular  waveguide  system 
then  extended  to  frequency-domain  and  time-domain  focus  arch  free-space  systems. 

Most  free-space  systems  consist  of  transmit  and  receive  spot-focusing  horn  lens 
antennas  and  a  network  analyzer  (NWA).  A  six  step  procedure  is  used  to  extract  the 
material  characteristics  using  the  NWA.  The  data  is  measured  (1)  using  the  NWA 
then  it  is  frequency  windowed  (2)  to  knock  down  the  sidelobes  that  appear  in  the 
time-domain  after  the  data  is  transformed.  Once  transformed  (3)  to  the  time-domain 
the  data  is  time  gated  (4)  to  eliminate  multiple  reflections  between  the  antennas  and 
sample.  The  data  is  transformed  (5)  back  into  the  frequency-domain  for  extraction 
(6)  of  permittivity  and  permeability  using  a  Nicolson-Ross-Weir  (NRW)  method  and 
two  transmission  method.  Windowing  and  transforming  the  data,  steps  two  and 
three,  decreases  accuracy  at  the  band  edges.  The  use  of  a  digital  oscilloscope  with  a 
time-domain  reflectometer  (TDR)  module  should  increase  accuracy  by  removing  the 
windowing  and  transforming  operations.  Thus,  a  secondary  goal  of  this  thesis  is  to 
analyze  the  effect  on  material  characterization  accuracy  using  an  oscilloscope  instead 
of  a  NWA  since  one  window  and  one  transform  operation  is  removed. 

Experimental  results  of  the  permittivity  and  permeability  values  of  a  mag¬ 
netic  radar  absorbing  material  (MRAM)  sample  for  the  three  measuring  devices  are 
presented.  The  two  transmission  method  is  compared  to  the  NRW  method  in  the 
waveguide  and  frequency-domain  focus  arch  system  to  validate  the  two  transmission 
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method.  A  differential  uncertainty  analysis  is  conducted  to  show  the  improved  accu¬ 


racy  of  the  two  transmission  method  over  a  reflection  measurement  based  parameter 
extraction  technique  such  as  the  NRW  method.  A  comparison  of  the  time-domain 
and  frequency-domain  focused  arch  results  is  also  conducted. 
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Investigation  of  Frequency-Domain 
and  Time-Domain  Free-Space 
Material  Measurements 

I.  Introduction 

Material  characterization  is  the  process  of  determining  the  relative  constitutive 
parameters  permittivity  and  permeability  (er,/j,r)  of  a  material.  The  permittivity  of 
a  material  determines  its  susceptibility  to  being  electrically  polarized  when  in  the 
presence  of  an  electric  field.  Likewise,  the  permeability  of  a  material  determines  its 
susceptibility  to  being  magnetically  polarized  when  in  the  presence  of  an  magnetic 
field.  Generally,  permittivity  and  permeability  are  complex  quantities,  where  the  real 
parts  are  related  to  energy  storage  and  the  imaginary  parts  are  related  to  energy  loss. 

Radar  absorbing  materials  have  many  uses  such  as  clutter  suppression  at  air¬ 
ports,  electromagnetic  shielding  used  in  electromagnetic  compatibility  (EMC)  appli¬ 
cations,  as  well  as  in  compact  ranges  and  anechoic  chambers.  The  Air  Force  uses 
radar  absorbing  material  (RAM)  on  many  of  its  low  observable  aircraft.  The  increas¬ 
ing  application  of  RAM  requires  greater  need  for  accurate  constitutive  parameter 
measurements. 

This  thesis  presents  a  new  transmission  only  material  measurement  technique 
and  differential  uncertainty  analysis.  The  technique  is  used  to  extract  er  and  /zr 
of  a  sample  and  is  first  demonstrated  in  a  waveguide  system.  Then  extended  to  a 
frequency-domain  and  time-domain  free-space  focus  arch  system  aimed  at  increasing 
measurement  accuracy  and  decreasing  processing  time  and  resources. 

1.1  Problem  Statement 

In  general  most  material  parameters  are  extracted  using  a  reflection  and  trans¬ 
mission  measurement.  The  use  of  a  reflection  measurement  can  introduce  large  error 
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due  to  the  position  dependence  of  the  reflection  measurement.  One  method  to  over¬ 
come  this  error  is  the  use  of  two  transmission  measurements  of  the  same  material  with 
different  thicknesses.  However,  an  argument  can  be  made  that  although  these  two 
pieces  of  material  with  differing  thicknesses  come  from  the  same  stock  they  are  indeed 
different  due  to  sample  fabrication  and  mixing  inconsistencies.  This  thesis  offers  a 
material  measurement  technique  that  uses  one  thickness  of  material  and  backs  it  by 
a  known  sample  thus  allowing  for  accurate  measurements  of  a  single  piece  of  material 
using  only  transmission  measurements. 

A  material  characterization  process  for  a  waveguide  system  consists  of  several 
basic  steps.  A  Thru-Reflect-Line  (TRL)  calibration  is  used  to  calibrate  the  waveguide 
system.  The  sample  must  be  accurately  machined  to  £11  the  waveguide  holder  to 
ensure  that  only  the  dominate  waveguide  mode  is  present.  A  NWA  is  used  to  measure 
the  reflected  and  transmitted  scattering  parameters  (S-parameters)  of  the  material. 
Numerical  algorithms  such  as  the  NRW  method  [1,  6,  8,  10,  11]  are  then  used  to 
extract  the  complex  permittivity  and  permeability  from  these  S-parameters.  The 
NRW  technique  requires  a  transmission  and  reflection  measurement  and  will  be  used 
to  show  the  error  caused  by  using  a  reflection  measurement.  If  the  NRW  algorithm 
cannot  be  used  an  iterative  method  such  as  a  Newton  Two-Dimensional  (2D)  root 
search  is  applied.  The  use  of  a  reflection  measurement  can  lead  to  large  error  due 
to  the  path  length  difference  that  can  occur  when  making  a  reflection  measurement 
as  shown  in  Figure  1.1.  A  position  correction  factor  and  the  use  of  both  forward 
and  reverse  measurements  can  be  used  to  account  for  this  error.  It  will  be  shown 
that  the  use  of  transmission  measurements  which  are  independent  of  sample  position 
eliminates  the  need  for  a  correction  factor  and  leads  to  more  accurate  results. 

The  material  characterization  process  for  a  focus  arch  free-space  system  is  some¬ 
what  different  from  the  waveguide  system.  The  focus  arch  system  consists  of  two 
wide  band  antennas,  two  dielectric  lenses,  a  sample  holder,  and  a  NWA.  The  sample 
is  placed  in  the  holder  between  two  dielectric  focusing  lenses  which  are  themselves 
between  two  wide  band  antennas.  The  NWA  is  then  connected  to  the  antennas  and 
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Figure  1.1:  Focus  arch  system  depicting  the  two  types  of  re¬ 
flection  measurement  error  that  can  occur.  One  error  source 
is  caused  by  the  sample  not  being  normal  to  the  incident  wave 
causing  a  increase  in  path  length.  The  second  error  source  is 
caused  by  the  sample  not  being  located  at  the  Z  —  0  reference 
plane. 
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the  S-parameters  of  the  material  are  measured.  The  lenses  serve  to  collimate  (equal 
phase  front)  and  focus  the  beam  on  the  material.  As  with  the  waveguide  system 
the  NRW  method  and  two  transmission  method  are  used  to  extract  the  er  and  pr 
from  the  measured  S-parameters.  Prior  to  applying  an  extraction  algorithm  a  cali¬ 
bration  process  must  be  conducted.  The  experimental  data  is  frequency  windowed 
to  knock  down  the  sidelobes  that  occur  when  transformed  to  the  time-domain.  This 
frequency  windowing  step  causes  a  loss  of  accurate  data  at  the  band  edges.  Once 
in  the  time-domain  the  data  is  time  gated  to  remove  unwanted  sample  antenna  in¬ 
teractions.  Finally,  the  data  is  transformed  back  into  the  frequency- domain  and  a 
numerical  algorithm  is  applied. 

The  material  characterization  process  for  a  focus  arch  free-space  system  can  be 
further  simplified  by  using  a  time-domain  system.  Replacing  the  NWA  with  a  digital 
oscilloscope  having  a  Time-Domain  Transmission  (TDT)  module  allows  for  direct 
time-domain  measurement  of  the  sample.  Eliminating  the  need  to  frequency  window 
and  transform  to  the  time-domain  improving  band  edge  data.  The  direct  time-domain 
measurements  are  then  gated  as  above  and  transformed  to  the  frequency-domain  for 
parameter  extraction. 

The  goal  of  this  thesis  is  to  first  demonstrate  the  two  transmission  technique 
in  a  rectangular  waveguide  system  eliminating  error  associated  with  the  reflection 
measurement.  The  transmission  method  is  then  extended  to  a  free-space  frequency- 
domain  focus  arch  system.  Lastly,  the  technique  is  further  extended  to  direct  time- 
domain  free-space  measurements  which  will  increase  band  edge  accuracy  and  reduce 
data  processing  time. 

1 . 2  Scope 

Many  other  measurement  devices  can  be  used  for  material  characterization  in¬ 
cluding,  but  not  limited  to,  circular  waveguide  systems,  stripline  systems,  and  coaxial 
waveguide  systems.  This  thesis  will  explore  S-band  (2.6-3.95  Gigahertz  (GHz))  and 
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X-band  (8.2-12.4  GHz)  rectangular  waveguide  applications  as  well  as  both  frequency- 
domain  and  time-domain  focus  arch  free-space  (2-18  GHz)  applications. 

1 . 3  Organization 

Chapter  2  provides  a  review  of  S-parameters  and  wave  transmission  matrices 
(A-parameters)  along  with  the  NRW  algorithm  for  a  waveguide  and  free-space  system. 
Chapter  3  presents  a  two  transmission  method  for  extracting  constitutive  parameters 
for  a  waveguide  system  and  free-space  system.  Chapter  4  presents  results  for  S- 
band  and  X-band  waveguide  transmission  measurements  and  an  analysis  of  frequency- 
domain  and  time-domain  free-space  results  is  also  included.  Finally,  conclusions  and 
recommendations  for  future  research  are  given  in  Chapter  5. 
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II.  Background 

This  chapter  provides  the  necessary  background  for  understanding  the  measurement 
methods  used  in  this  thesis. 


2.1  A-Parameter  Single  Layer  System 

A-parameters  are  a  simple  method  to  solve  single  and  multi-layered  systems. 
The  development  describes  the  relationship  between  incident  and  reflected  wave  am¬ 
plitudes  at  the  input  and  output  terminal  plane.  The  A-parameter  description  of  a 
single  layer  system  is  examined  first  then  generalized  to  a  N-layer  system  in  the  fol¬ 
lowing  section.  An  expanded  view  of  a  single-layered  environment  is  shown  in  Figure 
2.1.  If  Ri,  T12  and  R2,  T2\  are  the  respective  interfacial  reflection  and  transmission 
coefficients  due  to  waves  c\  and  b'2  then  the  following  relations  exist. 


c2  —  T\2C\  +  i?2&2 

b\  =  RiCi+T2ib2  (2.1) 


where  the  reflection  and  transmission  equations  are 


T\  —  TV2  —  1  +  R1 


T2 

R\ 


T21  —  1  +  R2 


-R'2 


Z2  +  Z\ 


Substituting  (2.2)  into  (2.1)  leads  to  the  matrix  expression 


Cl 

1 

1 

Ri 

C2 

. bl . 

“  Tx 

.  Rl 

1 

.  b'2  _ 

(2.2) 


(2.3) 


Equation  (2.3)  describes  the  relationship  between  the  incident  and  reflected  traveling 
waves  (ci,&i)  and  ( c2,&2 )  immediately  to  the  left  and  right  of  the  sample  interface, 
respectively.  It  is  necessary  to  show  how  the  waves  (c2,  b'2)  are  related  to  waves  (c2,  b2) 
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A  h 


Figure  2.1:  Single  layer  system  showing  the 

relationship  between  the  reflected  and  transmit¬ 
ted  coefficients  and  the  interfacial  reflection  and 
transmission  coefficients. 

located  a  distance  Z2  from  the  interface  as  shown  in  Figure  (2.1).  The  sample  region 
is  assumed  to  be  linear,  homogeneous,  and  isotropic  (simple  media)  therefore,  the 
following  relation  exists  between  the  waves. 


c2  =  c2e72Z2 
62  =  b2e~l2h 

This  relationship  can  be  written  in  matrix  form 


/ 

C2 

g72^2 

0 

C2 

b'2 

0 

g— 72^2 

b2 

(2.4) 


(2.5) 
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The  A-parameter  relationship  between  waves  (ci,fei)  and  (c2,b2)  is  obtained  by  sub¬ 
stituting  (2.5)  into  (2.3). 


Cl 

1 

g72^2 

Rie~^h 

C2 

An 

Ai2 

C2 

h 

“  Tx 

Riel2h 

g— 72^2 

b2 

A2i 

a22 

b2 

The  relationship  between  A-parameters  and  S-parameters  is  found  using  the  below 
equation 

bi  =  SuCi  +  S'1262  (2-7) 

c2  =  A21C1  +  S22b2  (2.8) 

Solving  (2.8)  for  the  incident  wave,  C\  gives  the  following 

Ci  =  yj~c2  -  |^62  (2.9) 

021  *J21 

where  the  coefficients  of  c2  and  b2  are  An  and  A12,  respectively.  Substituting  (2.9) 
into  (2.7)  gives  the  equation  for  the  reflected  wave  shown  below 

bi  =  An  ( ~  -^—b^  +  S\2b2  (2.10) 

\021  021  / 

Rearranging  and  grouping  terms  leads  to  the  below  equation  for  the  reflected  wave 


h 


11 


c2  + 


21 


An  S22  —  A12A21 


A 


21 


(2.11) 


where  the  coefficients  of  C2  and  b2  are  A2i  and  A22,  respectively.  The  relationship 
between  S-parameters  and  A-parameters  are  shown  below  in  matrix  form. 


An 

1 

CN 

1 

1 - 

(N 

1 

T— 1 

A21 

A22 

“  A^ 

ry  Q  Q  Q  Q 

on  021012  —  011022 
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(2.12) 


(2.13) 


1 

£o 

1 

<M 

to 

I 

1  1—1 

A  21  A11A22  —  A21A12 

_  S21 

i 

(N 

to 

'  An 

1  — a12 

Applying  the  above  relationship  to  (2.6)  leads  to  the  theoretical  forward  S-parameters 
of  a  single-layered  material 


nthy  _  A(1  —  P 2) 

11  “  1  —  R2P2 

nthy  _  -P(l  ~  R 2) 

21  “  1  —  R2P2 


(2.14) 

(2.15) 


where  R  is  the  reflection  coefficient  and  P  is  the  one-way  phase  delay  through  the 
material  given  by  [6],  [3]. 


P 

R 


Z2  +  Z  i 


(2.16) 


2.2  A-Parameter  Multi-Layered  System 

Figure  2.2  shows  a  multi-layered  system  in  free-space,  comprised  of  N  —  1  layers 
and  N  interfaces.  The  q,  bi  terms  are  the  incident  and  reflected  complex  wave  ampli¬ 
tudes  to  the  left  of  the  ith  interface.  The  length,  complex  permittivity,  and  complex 
permeability  of  the  ith  layer  are  denoted  by  k,  e+  and  /X;,  respectively.  The  interfacial 
reflection  and  transmission  coefficients  of  the  ith  interface  are  denoted  by  R{  and  T). 
From  the  single  layer  analysis  above  an  overall  A-parameter  description  of  a  cascaded 
N-la.yer  system  can  be  obtained  by  generalizing  (2.6) 


Cl 

7 

JV+l  , 

=  n+ 

bi 

•  1 
i=l 

Rie  lili 

o  ^Yi^i 


CN+1 

bN+i 


Asys  Asys 

^11  ^12 

Asys  Asys 

^21  ^22 


Rielik  e 

solving  for  the  system  A-parameters  leads  to  the  following  result 


Cn+ 1 

bN+i 

(2.17) 


A^s 

Asys 

JV+l 

=  n+ 

^TYi^-i 

Ri,e  lili 

JV+l 

=  n 

Ai 

Ai 

^12 

Asys 

A21 

A  SVS 
^22 

xx  Ti 

1=1 

Rielili 

g  'Yi^'i 

i= 1 

Ai 

^21 

Ai 

^22 

(2.18) 
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h 


k 


4r-i 


'i 


< - » 

c2  - ► 

4 - » 

bi  4 - 

b2  ◄ - 

Ti 

t2 

X" 

N  < 

N 

V> 

v 

(®o^o) 

(s2.^2) 

Figure  2.2:  A  multi-layered  system  in  free-space,  comprised  of  N  —  1  layers 
and  N  interfaces.  The  system  shows  the  relationship  between  the  reflected 
and  transmitted  coefficients  and  the  interfacial  reflection  and  transmission 
coefficients. 

where 


d  _  Aj—i 

‘  Zi  +  Zi_  i 

T%  —  1  +  Ri 
Pi  =  e~lik 


(2.19) 


As  an  example  consider  the  two-layered  sample  of  length  lsys  immersed  in  free-space 
shown  in  Figure  2.3.  There  are  three  interfaces  that  result  in  the  following  special¬ 
ization  of  (2.17) 


Cl 

3  1 

=lU 

elili  Rie  lik 

h 

Ti 

1=1 

1 

1 

03 

03 

c? 

a  svs  a  sys 

/ill  ^12 

A  sys  a  sys 

^21  ^22 


C4 

&4 


(2.20) 


^3  =  0  and  letting  P?;  =  e  lik  the  wave  matrix  for  the  system  is  therefore, 


Asys 


1 

T\T2T3 


pr1  RiPi 
RiPr1  Pi 


p~  1 
1  2 

R2P21 


R2P2 

P2 


(2.21) 
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V>s 


Figure  2.3:  Two-layered  sample  of  length  lsys  immersed  in  free-space  showing 
the  relationship  between  the  reflected  and  transmitted  coefficients  and  the 
interfacial  reflection  and  transmission  coefficients. 


where  ZQ 


representing  the  air  layer  [6] . 


R 


Z\  —  Z0 
1  “  ZX  +  Z0 
T\  —  1  +  R\ 


R2 


z  2  —  Z\ 
Z2  +  Z\ 


T2 

r3  = 
t3 


- 1  +  r2 
Zo  —  z2 
Zo  +  Z2 
=  1  +  R-.i 


(2.22) 
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2.3  NRW  Algorithm 

Material  measurements  compare  the  theoretical  S-parameters  (S^y,S^[y)  to  ex¬ 
perimentally  measured  sample  S-parameters  (Sl^p,S2iP) 

rithy.Ml /  \  riexp.Ml  /  \  n 

*321  (^,e,/i)-^2i  M  =  0 

S^’M2(uj,e,g)- SHp’M2(uj)  =  0  (2.23) 


Nicolson,  Ross,  and  Weir  combined  (2.9)  and  (2.10)  and  derived  explicit  formulas  for  a 
material’s  constitutive  parameters  er  and  gr.  Similarly  to  the  theoretical  S-parameter 
values  the  experimental  S-parameters  can  be  written  as 


^exp 


R(1  —  P2) 
1  -  R2P2 


(2.24) 


nexp 

°21 


P(  1  -  R2) 
1  -  R?P2 


(2.25) 


The  NRW  formulas  are  derived  by  first  solving  (2.24)  for  P 2  and  (2.25)  for  P 


P2 


R  -  SHP 
R{  1  -  RSHP) 


p  Se2 f(l  -  R2P2) 

1  -  R2 

Substituting  (2.26)  into  the  P 2  term  of  (2.27)  leads  to 


P 


Qexp 

*^21 


1  -  RSHP 


=*►  P 2 


Qexp 
^21 


RSHP 


(2.26) 


(2.27) 


(2.28) 


After  equating  (2.26)  and  (2.28) 

R  ~  Se2lp  _  f  S™p  \ 2  _  R  -  SHP  _  ( Se2lpf 
R{  1  -  RSHP)  ~\1-  RSI lp )  ^  R  “  1  -  RSHP 
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Simplifying  and  grouping  terms  leads  to  the  following  quadratic  equation 


R2  -  2  KR  +  1 


(2.30) 


where 


K 


(S™p)2  -  ( SHP)2  +  1 


(2.31) 


2SHP 


whose  solution  is 


R  =  K  ±  VK'2  -  1 


(2.32) 


The  correct  choice  of  positive  or  negative  sign  in  (2.32)  is  made  by  requiring  \R\  <  1. 
The  thickness  of  the  material  is  chosen  such  that  S^p  in  (2.31)  does  not  equal  zero. 
The  phase  delay  through  the  material  P  is  given  by  the  following 


The  below  equations  are  used  to  find  the  constitutive  parameters  (er,/xr)  of  a 
sample  of  length  l  in  a  rectangular  waveguide  operating  in  the  fundamental  transverse 
electric  (TE)  mode 


(2.34) 


(2.35) 


where 


(2.36) 
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The  constitutive  parameters  for  a  transverse  electromagnetic  (TEM)  system  can  be 
easily  found  by  letting  kc  —  0  reducing  the  above  waveguide  method  to  a  TEM 
method  [6],  [2]. 

2.f  Discrete  Fourier  Transform 

There  is  a  large  amount  of  data  processing  for  the  calibration  of  the  free-space 
system  including  transformations  to  and  from  the  frequency-domain  and  time-domain. 
Therefore,  the  discrete  fourier  transform  is  outlined  below  and  interested  readers  can 
obtain  a  further  review  of  Discrete  Fourier  Transforms  (DFT)  in,  ’Signal  Processing 
and  Linear  Systems,’  by  B.  P.  Lathi  [9]. 

The  following  equations  define  the  discrete  DFT  with  Fn  the  direct  DFT  of  f/,, 
and  fk  the  inverse  discrete  Fourier  transform  (IDFT)  of  Fn. 


N0-l 


k= 0 

Nr. —  1 


(2.37) 


The  notation  fk  <^>  Fn  denotes  a  DFT  pair.  Fn  is  the  nth  sample  of  F(u)  and  fk 
is  To/Nq  times  the  kth  sample  of  f{t).  If  the  sample  values  of  f(t )  are  known  then 


the  sample  values  of  F(u)  can  be  computed  and  vice  versa  using  the  DFT.  fk  is  a 
function  of  k  (k  —  0, 1,2,  ...No  —  1)  rather  than  that  of  t  and  Fn  is  a  function  of  n 
(n  =  0, 1,  2,  ...Nq  —  1)  rather  than  u.  A  sampled  signal  f(t)  can  be  expressed  as 


Nq-1 


m  =  V  f(kT)6(t  -  kT) 


(2.38) 


5(t-  kT)  e~iku* 


(2.39) 


(2.40) 
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The  Fourier  transform  of  (2.38)  yields 


JVo-l 

F(u)  =  E  f(kT)e~jkut 

k= 0 


(2.41) 


However,  over  the  interval  |cn|  <  y,  F(lu)  of  f(t )  is  assuming  negligible  aliasing 


N0-l 

F(u)  =  TF(uj)  =  T  E  f{kT)e~jk“T  (2.42) 

k= 0 


and 


iVo— 1 

Fn  =  F(nuj 0)  =  T  E  f(kT)e~jk™°T 

k= 0 


(2.43) 


Observing  the  fact  that  No,  the  number  of  samples  of  a  signal  in  one  period  To  is 
identical  to  Nq. 


No 

T0 

No 


To 

T 

1 

% 

To 

T 


=  N' 


Also,  substituting  Tf(kT)  =  /*,  yields  the  equation  for  the  DFT 


Wo— 1 

Fn=^2 

k= 0 


—jn^r-h 

e  No 


JmjFn 


To  find  the  inverse  relation  multiply  both  sides  of  (2.47)  by  eJ  No 


W0-l 


Em  jm-Tj^n 

FneJ  No  = 


n= 0 


Wo— 1 

E 

n=0 


W0  — 1 


E  f^k 


k= 0 


g  N0 


(2.44) 

(2.45) 

(2.46) 


(2.47) 


(2.48) 
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Interchanging  the  summation  on  the  right  hand  side 


No-l 

T.  Fne 


n= 0 


N0- 1 

y  fa 


k= 0 


~N0-1 

E 


ej{m-k)^n 


n= 0 


The  inner  sum  on  the  right  hand  side  is  zero  for  k  ^  m,  and  N0  when  k 
below 

N0-l 

£  =  N„gm,t 

n= 0 

Therefore,  N0fk  =  N0fm  which  yields  the  equation  for  the  IDFT 


JVo-l 


f  =  —  V  F  ejm^ 

Jm  Nn  ^  n 


n=0 


(2.49) 

m  shown 

(2.50) 


(2.51) 
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III.  Methodology 

This  chapter  presents  a  transmission  only  method  to  extract  the  constitutive  param¬ 
eters  of  a  test  sample  mounted  in  a  rectangular  waveguide  and  a  free-space  system. 
The  analysis  will  result  in  a  calculated  S-parameter  for  both  the  single-layer  sample, 
gthy,M i  an(j  £iie  two-layer  acrylic  backed  sample,  S%[y’M2.  The  theoretical  values  will 
be  compared  to  corresponding  single-layer  sample  measurements,  S^1’ M 1 ,  and  two- 
layer  acrylic  backed  sample  measurements,  S^[P’M2 .  Using  an  iterative  Newton  2D 
root  search  the  sample’s  material  parameters  will  be  those  theoretical  values  of  er,pr 
that  satisfy  the  below  equations  to  within  a  specified  accuracy  at  each  frequency: 

|d21  (a;,  er,  fir)  —  &2i  {co)\  <  accuracy 

\S^y,M2(u,er, /j,r)  —  S2iP,M2(uj)\  <  accuracy  (3.1) 

The  results  of  the  two  transmission  method  will  be  compared  to  the  NRW  method. 
It  will  be  shown  that  the  two  transmission  method  removes  path  length  dependent 
reflection  errors  associated  with  the  NRW  method. 

3. 1  Acrylic  Sample 

To  solve  for  both  er  and  /ir  it  is  necessary  to  have  two  independent  transmission 
measurements.  The  first  measurement,  ,  is  of  the  unknown  sample  only.  The 

second  experimental  measurement,  S^lP,M2  ^  is  the  sample  backed  by  a  piece  of  acrylic 
with  a  known  thickness  iacr  and  constitutive  parameters  e“cr  and  p,^cr.  Acrylic  is  a 
lossless  dielectric  having  a  relative  permeability,  yr  —  1  +  j 0.  This  fact  leads  to  a 
fundamental  frequency  phenomenon  when  the  phase  of  a  single-layer  material  given 
by  P  =  e~ikz°iacr  is  equal  to  a  multiple  of  a  half  wavelength.  At  the  fundamental 
frequency  when  m  =  1  there  is  no  phase  difference  and  very  little  magnitude  difference 
between  the  two  measurements.  The  phase  term  is  determined  by  the  waveguide 
propagation  constant  and  the  thickness  of  the  acrylic  sample 

9  =  kzJacr  (3.2) 
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where 


k::„  =  -  ki 


(3.3) 


The  free-space  propagation  constant  kQ  is  given  by 


kl  =  uj2e0h0 


(3.4) 


substituting  into  (3.3)  along  with  angular  frequency  u  =  2nf  results  in  the  following 


kz„  =  \J '2, -  k\ 


(3.5) 


The  above  equation  for  the  waveguide  propagation  constant  is  substituted  into  (3.2) 


V,  =  (y(2ir/)2£<,;i1>£“’>“r  -  fc?)  Cr  = 


(3.6) 


solving  for  /Q  results  in  the  following  expression  for  the  fundamental  waveguide  fre¬ 
quency 

I  <^  +  k2c 

fo=\  A  2  acr  acr  (3.7) 

V  47T 2£oldo£r  fj-r 

Setting  kc,  the  waveguide  cutoff  frequency,  equal  to  zero  and  applying  the  definition 
of  the  speed  of  light  in  free-space 


c  = 


\J  ^0^0 


(3.7)  simplifies  to  the  fundamental  free-space  frequency  given  below 


(3.8) 


/o  — 


2 


(3.9) 


At  the  fundamental  frequency  the  2D  root  search  will  not  converge  to  a  solution.  To 
ensure  the  two  measurements  have  sufficient  differences  in  phase  the  thickness  of  the 
acrylic  is  chosen  such  that  6  =  k0l  >  30°.  This  value  accounts  for  the  phase  accuracy 
of  the  NWA  which  is  roughly  15°  at  18  GHz  resulting  in  a  safety  margin  ensuring 
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convergence  to  a  solution.  A  maximum  phase  value  is  also  chosen  in  case  a  reflection 
measurement  is  taken  in  the  future.  In  the  case  of  a  reflection  measurement  the  phase 
is  a  two-way  value  and  a  180°  phase  difference  will  not  result  in  a  solution.  Therefore, 
the  maximum  phase  difference  is  chosen  to  be  6  =  <  150°  this  also  includes  a 

safety  margin.  This  phase  requirement  cannot  be  satisfied  over  the  entire  bandwidth 
of  the  focus  arch  system  (2-18  GHz)  using  a  single  piece  of  acrylic.  Therefore,  the 
band  is  broken  into  two  smaller  bandwidths  ranging  from  2-6  GHz  and  6-18  GHz 
requiring  two  pieces  of  acrylic.  The  thickness  is  required  to  be  |  at  mid-band  to 
ensure  optimum  incident  power  across  the  entire  band.  This  requirement  also  helps 
ensure  that  the  phase  across  the  band  does  not  approach  0°  or  180°.  The  wavelength 
A  in  a  rectangular  waveguide  is  given  by  the  following  equation 


A  = 


(3.10) 


Substituting  (3.5)  from  above  into  (3.10)  results  in  the  below  equation  for  wavelength 


A 


2tt 

/  I"-../'..-)' 7'C"  -  /.') 


(3.11) 


Requiring  £acr  =  |  and  solving  (3.11)  for  the  acrylic  thickness  at  the  mid  band 
frequency  fmid  gives  the  following 


yacr 

fmid 


2tt 


4v/(2vr/mi(i)2£0/i0£“c>rcr 


k 


(3.12) 


As  before  setting  kc  =  0  and  applying  the  definition  of  the  speed  of  light  in  free-space 
simplifies  (3.12)  to  give  the  acrylic  thickness  at  the  mid  band  frequency  for  free-space 
shown  below 


^ acr 
fmid 


4  fmids/£?Crl1 


acr 

r 


(3.13) 


The  acrylic  thickness  required  at  fmid  =  4  GHz  and  fmid  =  12  GHz  is  .453”  and 
.151”,  respectively.  Commercially  available  sizes  of  .5”  and  .125”  are  chosen  to  reduce 
cost.  The  minimum  and  maximum  phase  across  the  band  is  verified  using  (3.2). 
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Substituting  (3.5)  into  (3.2)  leads  to  the  equation  for  the  waveguide  phase  at  the 
minimum  and  maximum  frequencies  fmin,max  of  the  band 


9  = 


\j  (2vr/, 


min, max , 


ic-  ii  c-acri,acr 

^oH'oC'r  r^r 


(3.14) 


Setting  kc  =  0  and  applying  the  definition  of  the  speed  of  light  in  free-space  simplifies 
(3.14)  to  give  the  minimum  and  maximum  free-space  phase 


9 


f ‘min,  max  £  £rfir 


C 


(3.15) 


For  the  .5”  acrylic  sample  the  phase  across  the  band  ranges  from  45.6°  <  9  <  137.1° 
and  for  the  .125”  sample  the  phase  ranges  from  34.3°  <  9  <  102.8°.  The  phase 
difference  for  both  waveguide  systems  can  be  easily  satisfied  using  the  above  acrylic 
thicknesses. 


3.2  Waveguide 

Both  S-band  (2.6-3.95  GHz)  and  X-band  (8.2-12.4  GHz)  rectangular  waveguides 
are  used  to  demonstrate  the  transmission  technique.  Figure  3.1(a)  shows  the  S-band 
waveguide  including  the  measured  samples  and  the  sample  holder.  The  cutoff  wave 
number  for  the  S-band  waveguide  is  given  by 

kc=-  (3.16) 

a 

where  a  is  the  longest  dimension  of  the  waveguide,  for  S-band  a  =  .072136  m.  The 
X-band  waveguide,  measured  samples,  and  sample  holder  are  shown  in  Figure  3.1(b) 
whose  longest  dimension  is  a  =  .02286  m. 

3.2.1  Calibration.  Before  taking  any  measurements  it  is  necessary  to  cal¬ 
ibrate  the  waveguide  system  to  remove  any  system  errors  and  establish  the  phase 
reference  point.  This  Z  —  0  phase  reference  plane  is  the  location  of  the  front  face  of 
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Figure  3.1:  (a)  S-band  (2.6-3.95  GHz)  rectangular  waveguide  shown  with  MRAM 

sample,  acrylic  sample,  and  sample  holder. 

(b)  X-band  (8.2-12.4  GHZ)  rectangular  waveguide  shown  with  MRAM  sample,  acrylic 
sample,  and  sample  holder. 


the  sample.  The  sample  will  be  located  at  this  location  for  all  measurements  to  ensure 
the  phase  delay  term  is  calculated  correctly.  Figure  3.2  shows  the  three  calibration 
measurements  necessary  to  calibrate  the  waveguide  system.  The  HP8510  NWA  must 
be  setup  prior  to  calibration.  For  S-band  measurements  the  start  frequency  2.6  GHz 
and  stop  frequency  3.95  GHz  are  set  along  with  the  waveguide  cutoff  frequency  (re¬ 
ferred  to  as  waveguide  delay  in  the  NWA)  of  2.078  GHz.  Similarly,  for  the  X-band 
measurements  the  start  frequency  8.2  GHz  and  stop  frequency  12.4  GHz  are  set  along 
with  the  waveguide  cutoff  frequency  of  6.557  GHz.  The  following  settings  are  common 
between  both  waveguide  systems.  The  NWA  is  set  to  measure  201  frequency  points 
with  a  100  ms  sweep  time  taking  64  averages  at  each  frequency.  Once  the  NWA  pa¬ 
rameters  are  set  the  TRL  calibration  is  conducted  using  a  zero  delay  for  the  sample 
holder.  This  zero  delay  accommodates  the  holder  for  the  measurement  eliminating 
the  need  to  take  the  length  of  the  holder  into  account.  The  delay  for  the  S-band  and 
X-band  line  standards  are  -32.376  ps  and  -101.53  ps,  respectively.  The  first  calibra¬ 
tion  measurement  is  the  thru  measurement  with  the  sample  holder  attached.  The 
next  calibration  measurement  is  the  reflect,  made  by  inserting  a  short  as  depicted 
in  Figure  3.2.  The  final  calibration  measurement  is  the  line  measurement  made  by 
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z=o  z=o 


Figure  3.2:  Depiction  of  the  calibration  setup  used  for  both 

S-band  and  X-band  waveguide  TRL  calibration. 

connecting  the  waveguide  without  the  sample  holder.  After  omitting  the  isolation 
calibration  measurement  and  saving  the  calibration  coefficients,  it  is  good  practice  to 
measure  S\  i  to  ensure  it  has  a  magnitude  of  one  and  a  phase  of  180°. 


3.2.2  Waveguide  NRW  Method.  The  sample  is  cut  to  fill  the  entire  waveg¬ 
uide  holder  to  ensure  dominant  mode  operation.  Forward  reflection  S '“P’measl;  anc[ 
transmission  g^P’measl  S-parameters  are  measured  and  recorded.  Prior  to  parameter 
extraction  the  measured  S-parameters  must  be  phase  shifted  to  account  for  the  thick¬ 
ness  of  the  sample  because  it  is  absent  in  the  thru  measurement  as  shown  in  Figure 
3.3(a).  The  S-parameters  are  shifted  using  the  below  equations 


gexp,M  1 

_  gexp,measl 

Qexp,M  2 
*^21 

_  gexp,meas2  ^-jkZoes 

(3.17) 
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Figure  3.3:  (a)  Depiction  of  phase  relationship  between  measured  S-parameters  and 

measured  sample  S-parameters  for  a  single-layer  system. 

(b)  Depiction  of  phase  relationship  between  measured  S-parameters  and  measured 
sample  S-parameters  for  a  two-layer  system. 

where  kZo  is  the  propagation  constant  and  is  is  the  sample  thickness.  The  above 
experimental  S-parameters  are  used  in  the  NRW  algorithm  outlined  in  Chapter  2 
to  extract  er  and  fir.  The  extracted  complex  relative  permittivity  ereai  +  jeimag  and 
complex  relative  permeability  /irea;  +  jfiimag  of  a  MRAM  sample  in  a  S-band  and 
X-band  rectangular  waveguide  are  presented  in  chapter  four. 

3.2.3  Waveguide  Two  Transmission  Method.  The  two  transmission  method 
uses  the  above  single-layer  sample  measurement,  S '^p’measl)  and  a  two-layer  acrylic 
backed  sample  measurement,  S% ^p’meas2_  Measurement  one  is  phase  shifted  as  shown 
in  Figure  3.3(a)  while  measurement  two  is  phase  shifted  as  shown  in  Figure  3.3(b). 
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The  S-parameters  are  shifted  using  the  below  equations 


Qexp,M  1 
°21 

Qexp,M  2 
*~^21 


gexp, measly  -j  kZo  t3 
gexp,meas2  g  -  j  kZo  ( ls  +£acr ) 


(3.18) 


where  kZo  is  the  propagation  constant,  £s  is  the  sample  thickness,  and  £acr  is  the  acrylic 
thickness.  These  experimental  values  are  compared  to  calculated  theoretical  values 
using  a  2D  Newton  root  search  solving  the  below  equations 

|o21  (u>,  e,  /i)  —  o21  (cn)|  <  accuracy 

\St2^y,M2(uji  e,  n)  —  S^[P,M‘2  {u)\  <  accuracy  (3.19) 


An  initial  guess  is  made  for  the  unknown  sample  er  and  /ir  and  the  root  search 
iteratively  approaches  a  solution  that  satisfies  (3.19)  to  within  a  specified  accuracy. 
The  theoretical  transmission  coefficient  for  the  single-layer  measurement  S'2p/,i'n  is 
given  by  2.10  and  shown  below 


Qthy,M  1 
*^21 


P(1  -R2) 
1  -  R2P2 


(3.20) 
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The  theoretical  transmission  coefficient  for  the  two-layer  measurement  is  found  by 
solving  (2.21)  for  An  and  using  (2.8)  to  solve  for  S2i 


nthy,M  2 
^2l' 


_ TlT2T3 _ 

(Pf^-1  +  R.RzP.P,-1)  +  (Pf1P2P2  +  R1P1P2)R3 


(3.22) 


where  Ri,  T),  and  P*  are  the  reflection  coefficient,  transmission  coefficient,  and  one¬ 
way  phase  delay  for  the  ith  interface.  The  theoretical  transmission  coefficient  for  a 
rectangular  waveguide  can  be  further  simplified  to  the  following  [7] 


Qthy,M  2 
D21 


Y± 


8Z0ZUnk 


e7 acrtacr  [Z acr  +  Z Q)  Y+  +  eT^lacr  (Z Q  -  Z acr)  Y~ 

2Zunj^Z acr  cos  kZUnk^s  i  2jZunk  sin  k  zunk^s  Y  2jZ0Zacr  sin  A 'zunk^s 


(3.23) 


it  2Z0Zunk  COS  k zunk f  .s 


where  Zunk  is  the  waveguide  impedance  of  the  unknown  layer,  Zacr  is  the  waveguide 
impedance  of  the  acrylic  layer,  and  kzunk  is  the  propagation  constant  in  the  unknown 
layer.  The  extracted  complex  relative  permittivity  ereai  +  jeimag  and  complex  relative 
permeability  fireai  +  jUimag  of  a  MRAM  sample  in  a  S-band  and  X-band  rectangular 
waveguide  are  presented  in  chapter  four. 


3.3  Frequency- Domain  Focus  Arch 

The  frequency-domain  focus  arch  system,  free-space  system,  consists  of  a  NWA 
connected  to  a  transmit  horn  antenna  through  which  a  TEM  wave  is  transmitted. 
The  wave  passes  through  a  system  of  dielectric  lenses  that  collimate  and  focus  the 
energy  onto  the  sample.  The  wave  then  passes  through  or  is  reflected  by  the  sample 
and  is  received  by  another  horn  antenna  as  shown  in  Figure  3.4.  Since  the  sample  size 
is  large  relative  to  wavelength  and  the  lens  system  focuses  the  beam  to  a  spot  on  the 
sample,  it  can  be  approximated  as  infinite  in  extent  and  edge  scattering  effects  can 
be  ignored.  Unlike  the  waveguide  procedure  which  uses  a  TRL  calibration  the  focus 
arch  system  requires  a  significant  amount  of  postprocessing  to  calibrate  the  system. 
The  calibration  procedure  consists  of  a  simple  response  calibration  and  range-gating 
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Figure  3.4:  General  Electric’s  focus  arch  free-space  range  with 
transmit  horn  on  the  left  and  receive  horn  on  the  right.  The 
sample  holder  is  not  shown  but  is  located  in  the  middle,  in  place 
of  the  calibration  cylinder  shown. 

(time-gating)  to  calibrate  the  system.  Once  calibrated  the  NRW  method  and  two- 
transmission  root  search  method  are  used  to  extract  er,/ir. 

3.3.1  Two-Transmission  Calibration  and  Extraction.  The  two  transmission 
method  requires  three  measurements  a  single-layer  sample  measurement  S '“P>measl;  a 
two-layer  acrylic  backed  sample  measurement  g^P’meas2^  anc[  an  empty  measurement 
gexp, empty '  ^iree  measurements  are  processed  the  same  but  for  illustrative  pur¬ 

poses  only  g^P’measl  data  is  shown.  Calibration  of  the  free-space  frequency-domain 
system  is  conducted  by  first  windowing  the  S '^p’measl  data  shown  in  Figure  3.5(a). 
The  windowing  is  necessary  due  to  the  finite  bandwidth  of  the  measurement.  When 
transformed  to  the  time-domain  the  dynamic  range  of  the  time-domain  measurement 
is  limited  by  hiding  low  level  responses  within  the  higher  level  sidelobes.  Frequency 
windowing  the  measured  data  serves  to  lower  the  sidelobes  that  occur  when  trans¬ 
formed  to  the  time-domain.  The  Kaiser-Bessel  window  with  (3  =  6  is  chosen  and 
reduces  the  sidelobes  50  Decibels  (dB)  relative  to  the  peak.  The  frequency-windowed 
data  is  shown  in  Figure  3.5(b).  The  windowing  function  causes  a  rounding  ef¬ 
fect  of  the  band  edges  reducing  valid  data  at  the  band  edges.  Once  windowed  and 
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(a)  (b) 

Figure  3.5:  (a)  Raw  measured  forward  transmission  measurement  S '^p’measl  prior 

to  calibration  and  normalization  of  a  single-layer  sample  in  the  focus  arch  system. 

(b)  Windowed  forward  transmission  measurement. 


transformed  to  the  time-domain  as  shown  in  Figure  3.6(a)  any  unwanted  interaction 
between  the  sample  and  lens  system  is  gated  out.  The  time  gate  is  applied  1  ns  around 
the  location  of  the  sample  to  ensure  lens  sample  interaction  is  removed  as  shown  in 
Figure  3.6(b).  After  time  gating  the  data  is  transformed  back  to  the  frequency-domain 
and  is  normalized  to  the  empty  measurement  and  phase  delayed  to  the  Z  =  0  plane 
shown  in  Figure  3.3.  The  equations  for  the  simple  response  calibration  and  phase 
delay  are  given  by 

Qexp,Ml 
^21 

Qexp,M  2 
*~>21 

where  k0  is  the  free-space  propagation  constant,  is  is  the  sample  thickness,  laCr  is 
the  acrylic  thickness,  and  S^p’ernpl u  is  the  forward  transmission  measurement  with  no 
sample  in  the  holder.  The  resulting  calibrated  data  shown  in  Figure  3.7  is  used  in  a 
2D  root  search  that  compare  experimental  values  to  the  calculated  theoretical  values 
using  (3.19).  The  extracted  complex  relative  permittivity  ereai  +  jtimag  and  complex 
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(3.24) 
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Figure  3.6:  (a)  Results  of  IDFT  on  windowed  forward  transmission  measurement, 

(b)  Gated  time-domain  forward  transmission  measurement. 

relative  permeability  nreai  +  j^imag  of  a  MRAM  sample  in  a  focus  arch  system  are 
presented  in  chapter  four. 


3.3.2  NRW  Calibration  and  Extraction.  Unlike  the  two  transmission  method 
the  NRW  extraction  method  requires  four  measurements  a  short  reflection  S'^p’s/lort, 
empty  transmission  o21  ,  sample  reflection  on  ,  and  sample  transmission 

^exp, sample  measuremenp  The  four  measurements  are  processed  as  stated  above  to 
remove  any  sample  lens  interactions.  A  simple  response  calibration  and  phase  delay 
is  applied  given  by 


SexP 


nexp 

^21 


gexp 


,  short 


e~jk  o  (la+l  acr') 


(3.25) 


The  above  calibrated  experimental  S-parameters  are  used  in  the  NRW  algorithm 
outlined  in  Chapter  2  to  extract  er,/ir.  The  extracted  complex  relative  permittivity 
Ceai  +  jCmag  and  complex  relative  permeability  /irea;  +  jinimag  of  a  MRAM  sample  in 
a  focus  arch  system  are  presented  and  compared  to  the  two  transmission  method  in 
chapter  four. 
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Figure  3.7:  Depiction  of  focus  arch  system  for  a  single-layer 

system. 

3.4  Direct  Time- Domain  Focus  Arch 

The  time-domain  focus  arch,  free-space  system,  consists  of  an  Agilent  86100B 
digital  oscilloscope  with  54754A  differential  time-domain  transmission  (TDT)  module 
connected  to  a  power  amplifier  that  is  connected  to  a  transmit  horn  antenna.  The 
TDT  module  within  the  Oscilloscope  creates  a  200  mV  step  with  a  40  ps  rise  time  as 
shown  in  Figure  3.9.  The  antenna  removes  any  DC  component  from  the  step  input 
and  transmits  all  frequencies  within  the  spectral  content  of  the  step  input.  Unlike  the 
frequency-domain  focus  arch  system  the  calibration  procedure  is  simplified  eliminating 
the  frequency  window  and  transformation  to  the  time-domain  steps.  The  use  of  direct 
time-domain  measurements  serves  to  increase  the  band  edge  data  by  not  requiring  a 
window  function.  Once  calibrated  the  two-transmission  root  search  method  is  used  to 
extract  er  and  /ir.  The  NRW  algorithm  is  not  used  to  extract  constitutive  parameters 
in  the  case  of  the  time-domain  free-space  system  due  to  the  difficulty  in  taking  a 
reflection  measurement.  Due  to  the  low  output  voltage  of  the  TDR/TDT  module  the 
antenna  acts  as  a  reactive  load  and  does  not  allow  enough  energy  to  propagate  and  a 
received  signal  cannot  be  measured.  With  the  use  of  a  power  amplifier  the  antenna  is 
able  to  propagate  a  measurable  signal.  However,  the  use  of  a  power  amplifier  causes 
difficulty  in  making  a  reflection  measurement.  A  power  divider  and  an  additional  cable 
is  needed  to  receive  a  reflection  signal.  Additionally,  a  time  delay  occurs  between  the 
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Power  Supply 


Figure  3.8:  Depiction  of  time-domain  focus  arch  system  using 
a  digital  oscilloscope  with  TDR/TDT  module  connected  to  a 
power  amplifier. 


reflection  and  transmission  measurement  due  to  the  power  divider  and  must  also  be 
taken  into  account  to  ensure  correct  phase  calculation.  A  reflection  measurement  can 
be  made  with  the  time-domain  system,  particularly  in  systems  that  do  not  require  an 
amplifier.  However,  time-domain  reflection  measurements  in  a  focus  arch  free-space 
system  are  difficult  and  are  omitted  in  this  thesis. 

3-4-1  Two-Transmission  Calibration  and  Extraction.  The  time-domain  two 
transmission  method  requires  three  measurements  a  single-layer  sample  measurement 
V  {t)2 1  ,  a  two-layer  acrylic  backed  sample  measurement  V  (rj2i  >  and  an 

empty  measurement  V(t) ^iJ>emPty _  All  three  measurements  are  processed  the  same 
but  for  illustrative  purposes  only  V(t) ®  ‘p,measl  data  is  shown.  Calibration  of  the  free- 
space  time-domain  system  is  conducted  by  first  time  gating  the  y {t)'^p,rneasl  measured 
data  shown  in  Figure  3.10(a).  It  should  be  noted  that  the  derivative  of  the  measured 
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Figure  3.9:  Digital  oscilloscope  TDT  module  output  showing 
the  200  mV  step  with  a  40  ps  rise  time. 

data  may  need  to  be  taken  prior  to  time  gating  when  a  DC  component  is  present. 
For  instance  when  using  a  stripline  device  the  derivative  of  the  step  response  must 
be  taken  to  remove  the  DC  component  allowing  for  easier  time-gating.  Due  to  the 
nature  of  the  focus  arch  system,  the  DC  component  is  removed  by  the  antenna  and 
the  derivative  is  not  necessary.  Time  gating  becomes  much  more  of  an  art  than  a 
science  in  this  case.  For  the  focus  arch  system  the  time  gating  is  applied  at  a  zero 
point  in  the  step  response.  If  the  gate  is  applied  at  a  non-zero  point  a  rounding  of 
the  data  at  the  band  edges  can  occur.  To  ensure  valid  data  across  the  entire  band  the 
time  gate  must  end  at  a  zero  crossing  as  shown  in  Figure  3.10(b).  Once  time  gated 
the  data  is  transformed  to  the  frequency-domain  and  a  simple  response  calibration 
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Figure  3.10:  (a)  Raw  measured  forward  voltage  measurement  V (t) ^p,measl  using 

oscilloscope  prior  to  calibration  of  a  single-layer  sample  in  time-domain  focus  arch 
system. 

(b)  Time  gated  forward  voltage  measurement  V (£)^yi’Tneasl  prior  to  calibration  of  a 
single-layer  sample  in  time-domain  focus  arch  system. 


and  phase  delay  is  applied  given  by 
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where  the  transmission  S-parameter  is 
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Figure  3.11:  Calibrated  experimental  S-parameters  used  in 

2D  Newton  root  search  to  extract  er  and  f. ir. 

resulting  in 


Qexp,Ml 

o2i  — 

\r ( .  ,\trans,meas  1 

17  C-M 

Y  an  s,  empty 

q  exp,  M2 
°21  — 

,\trans,meas2 

V  c-ikM„+taar) 

V  (^jj^rans->ernP^y 

(3.28) 


As  shown  in  (3.28)  the  incident  voltage  cancels  out  and  a  step  or  impulse  input  can 
be  used.  The  above  calibrated  experimental  S-parameters  shown  in  Figure  3.11  are 
used  in  the  2D  Newton  root  search  to  extract  er  and  /ir.  The  extracted  complex 
relative  permittivity  ereai  +  jeirnag  and  complex  relative  permeability  fireai  +  j^imag  of 
a  MRAM  sample  in  a  time-domain  focus  arch  system  are  presented  and  compared  to 
the  frequency-domain  focus  arch  results  in  chapter  four. 
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IV.  Results  and  Analysis 

This  chapter  presents  the  results  of  the  two  transmission  method  for  a  rectangular 
waveguide  system,  frequency-domain  focus  arch  measurement  system,  and  a  time- 
domain  focus  arch  measurement  system.  The  complex  permittivity  and  permeability 
of  a  40  mil  FGM-40  MRAM  sample  are  presented  along  with  a  differential  uncertainty 
analysis.  The  sample  is  cut  to  fill  the  entire  cross  sectional  area  of  the  waveguides 
and  the  free-space  sample  is  a  12  in  square.  The  two  transmission  method  can  be 
applied  to  non-magnetic  materials.  A  2D  root  search  is  not  necessary  for  non-magnetic 
materials  however,  the  two  transmission  method  can  be  used  and  the  results  for  a  125 
mil  Electric  Radar  Absorbing  Material (ERAM) sample  are  presented  in  Appendix  A. 

4-1  Waveguide  System 

4-1.1  Uncertainty  analysis.  Uncertainties  in  the  experimental  S-parameters 
will  lead  to  different  solutions  that  satisfy  (3.1).  There  are  many  sources  of  error 
that  contribute  to  the  uncertainty  in  the  waveguide  measurement  of  er  and  jir.  These 
include  alignment  errors  causing  higher  order  modes,  uncertainties  in  the  sample 
dimensions,  and  absorption  due  to  imperfectly  conducting  waveguide  walls  [4].  As¬ 
suming  the  walls  are  perfectly  conducting  and  the  sample  fully  fills  the  cross  section 
of  the  waveguide.  The  main  sources  of  uncertainty  for  the  NRW  method  are  the  un¬ 
certainty  of  the  sample  thickness  Sample  and  reflection  path  length  difference  x.  The 
main  sources  of  uncertainty  for  the  two  transmission  method  are  the  uncertainty  in 
the  sample  thickness  l sample,  acrylic  thickness  (-acrylic,  and  the  relative  permittivity  of 
the  acrylic  e“cr.  The  sample  thickness  is  varied  by  (sample  =  ±1%  and  the  path  length 
difference  is  varied  by  a  nominal  value  of  x  —  ±5  mils  for  the  NRW  method.  For 
the  transmission  method  the  sample  thickness,  acrylic  thickness,  and  acrylic  permit¬ 
tivity  are  varied  by  ±1%.  The  2D  root  search  and  NRW  method  are  run  with  these 
measurement  uncertainties  and  the  bounds  are  produced  as  shown  in  Figure  4.2. 

4-1.2  S-band  (2.6-3.95  GHz).  Figure  4.2  shows  the  complex  permittivity 
of  a  MRAM  sample  using  a  S-band  rectangular  waveguide  for  the  forward  NRW  and 
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two  transmission  method.  As  can  be  seen  from  the  figure  the  use  of  the  reflection 
measurement  in  the  NRW  method  causes  a  large  uncertainty.  The  two  transmission 
method  is  less  sensitive  to  the  position  of  the  material  and  therefore  has  a  much  smaller 
uncertainty  and  comparable  values  to  the  NRW  method.  Some  of  the  uncertainty  in 
the  measurements  may  be  due  to  the  low  frequencies  and  thin  sample.  The  complex 
permeability  of  the  MRAM  sample  is  shown  in  Figure  4.3.  The  results  for  the  S-band 
rectangular  waveguide  shows  the  two  transmission  method  is  as  good  or  better  than 
the  NRW  method  with  a  lower  opportunity  for  uncertainty. 

4- 1.3  X-band  (8.2-12-4  GHz).  The  complex  permittivity  of  the  MRAM  sam¬ 
ple  using  a  X-band  rectangular  waveguide  for  the  forward  NRW  and  two  transmission 
methods  are  shown  in  Figure  4.4.  Similar  to  the  S-band  measurements  the  X-band 
measurements  exhibit  a  large  uncertainty  when  using  the  NRW  method.  This  uncer¬ 
tainty  is  increased  in  the  X-band  waveguide  due  to  the  higher  frequencies.  The  two 
transmission  method  exhibits  a  much  smaller  uncertainty  and  comparable  values  to 
the  NRW  method.  The  complex  permeability  of  the  MRAM  sample  is  shown  in  Fig¬ 
ure  4.5.  The  results  for  the  X-band  rectangular  waveguide  shows  the  two  transmission 
method  is  more  accurate  measurement  method. 

4-2  Focus  Arch  System  (2-18  GHz) 

4-2.1  Uncertainty  analysis.  Some  of  the  errors  in  the  free-space  system  are 
the  collimation  and  focusing  of  the  beam,  sample  antenna  interactions,  and  back¬ 
ground  interference.  These  errors  are  predominantly  eliminated  through  the  simple 
response  and  range  gating  calibration.  The  sample  thickness  and  path  length  differ¬ 
ence  are  the  main  sources  of  uncertainty  for  the  NRW  method.  Sample  thickness  is 
varied  by  i sample  =  ±1%  and  the  path  length  difference  is  varied  by  a  nominal  value 
of  x  =  ±10  mils  for  the  NRW  method.  The  increase  in  the  reflection  path  length 
uncertainty  is  due  to  the  increase  in  path  length  and  sample  size.  The  path  length  for 
the  waveguide  system  is  less  than  a  foot  whereas  the  free-space  path  length  is  6  feet 
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for  the  free-space  system.  The  waveguide  system  offers  greater  accuracy  in  sample 
placement.  For  the  transmission  method  the  sample  thickness,  acrylic  thickness,  and 
acrylic  permittivity  are  varied  by  ±1%.  The  2D  root  search  and  NRW  method  are 
run  with  these  measurement  uncertainties  and  the  bounds  are  produced  as  shown  in 
Figure  4.6. 

4-2.2  Frequency-Domain.  Figure  4.6  shows  the  complex  permittivity  of  a 
MRAM  sample  in  a  focus  arch  system  from  2-6  GHz  using  the  forward  NRW  and 
two  transmission  methods.  As  can  be  seen  from  the  figure  the  use  of  the  reflection 
measurement  in  the  NRW  method  causes  a  slightly  larger  uncertainty  than  the  two 
transmission  measurements.  Both  extraction  methods  compare  well  with  the  values 
found  using  the  waveguide  system.  Note  the  band  edge  roll  off  that  occurs  at  2.5  and 
5.5  GHz.  The  complex  permeability  of  the  MRAM  sample  is  shown  in  Figure  4.7. 
The  results  for  the  frequency-domain  focus  arch  system  shows  the  two  transmission 
method  has  comparable  values  to  the  NRW  method  with  lower  uncertainty.  The 
complex  permittivity  of  a  MRAM  sample  in  a  focus  arch  system  from  6-18  GHz  using 
the  forward  NRW  and  two  transmission  methods  are  shown  in  Figure  4.8.  As  can 
be  seen  from  the  figure  the  use  of  the  reflection  measurement  in  the  NRW  method 
causes  a  much  larger  uncertainty  than  the  two  transmission  measurements  due  to  the 
higher  frequencies  (smaller  wavelengths).  Again  values  for  both  extraction  methods 
compare  well  with  the  values  found  using  the  waveguide  system.  Note  the  band 
edge  roll  off  that  occurs  at  6.5  GHz  and  17  GHz  due  to  windowing  the  frequency 
data.  The  complex  permeability  of  the  MRAM  sample  is  shown  in  Figure  4.9.  The 
results  at  higher  frequencies  for  the  frequency-domain  focus  arch  system  shows  the 
two  transmission  method  is  a  more  accurate  measurement  method  than  the  NRW 
algorithm. 

4-2.3  Time-Domain.  The  complex  permittivity  of  a  MRAM  sample  in  a 
frequency-domain  and  time-domain  focus  arch  system  from  2-6  GHz  using  the  two 
transmission  method  are  shown  in  Figure  4.10.  The  uncertainty  in  the  two  measure- 
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Figure  4.1:  Comparison  of  the  normalized  spectral  content 

contained  in  a  40  ps  step  input  signal  from  the  TDR/TDT  mod¬ 
ule  and  a  15  ps  step  input  signal  from  a  pulse  generator. 

ment  systems  are  calculated  the  same.  The  real  difference  between  these  two  methods 
is  the  increase  in  valuable  band  edge  data.  The  time-domain  results  hold  constant 
throughout  the  band  while  the  frequency-domain  results  diverge  at  the  edges.  The 
complex  permeability  of  the  MRAM  sample  is  shown  in  Figure  4.11.  The  time-domain 
results  in  the  lower  frequency  band  is  comparable  to  the  frequency- domain  results. 
Figure  4.12  shows  the  complex  permittivity  of  a  MRAM  sample  in  a  frequency-domain 
and  time-domain  focus  arch  system  from  6-18  GHz  using  the  two  transmission  method. 
The  complex  permeability  of  the  MRAM  sample  is  shown  in  Figure  4.13.  The  time- 
domain  measurements  increase  band  edge  data  as  expected.  However,  at  frequencies 
above  10  GHz  the  time-domain  results  diverge  from  the  frequency-domain  results. 
This  is  due  to  the  limited  spectral  content  of  the  40  ps  rise  time  step  input.  To  in¬ 
crease  accuracies  at  higher  frequencies  it  will  be  necessary  to  increase  the  rise  time  of 
the  step  input.  Figure  4.1  shows  the  normalized  spectral  content  of  the  time-domain 
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system  using  a  200  mV  step  input  compared  to  a  -5  V  15  ps  step  input  using  a  pulse 
generator.  Note  the  decrease  in  spectral  content  at  higher  frequencies  resulting  from 
the  slow  rise  time  of  the  step.  The  magnitude  of  the  spectral  content  at  10  GHz 
is  half  that  of  the  2  GHz  value.  Figure  4.1  also  shows  the  increased  spectral  con¬ 
tent  achieved  by  using  a  faster  rise  time  (fall  time)  using  a  pulse  generator.  Overall 
the  direct  time-domain  results  compare  well  with  the  frequency-domain  but  is  most 
accurate  at  lower  frequencies. 

4-2-4  Frequency- Domain  and  Time-Domain  Comparison.  Some  of  the  bene¬ 
fits  in  using  the  frequency-domain  system  include  a  stable  signal  source,  large  dynamic 
range,  and  lower  external  interference.  Benefits  of  the  time-domain  system  include 
easily  gating  multipath  interference,  broadband  data  is  obtained  in  one  measurement, 
and  measurements  are  very  rapid.  There  are  also  drawbacks  to  each  system.  The 
frequency-domain  system  requires  large  averaging,  frequency  windowing,  and  trans¬ 
formation  to  the  time-domain  is  required  for  gating.  The  time-domain  system  requires 
large  averaging  to  achieve  good  Signal-to-Noise  Ratio  (SNR),  fast  rise-time  is  required 
for  high  spectral  content,  and  it  is  prone  to  interference  from  external  sources  [5].  The 
direct  time-domain  system  requires  a  power  amplifier  with  a  fast  response  time  which 
can  be  quite  costly.  The  time-domain  system  also  requires  additional  hardware  to 
achieve  higher  frequency  measurements. 

4  •  3  Summary 

The  results  show  that  the  two  transmission  method  is  a  valid  measurement 
method  for  measuring  the  complex  permittivity  and  permeability  of  a  MRAM  sample. 
Furthermore,  the  two  transmission  method  has  a  lower  uncertainty  than  the  NRW 
method  especially  at  higher  frequencies.  Using  the  two  transmission  method  in  a 
direct  time-domain  focus  arch  system  can  further  increase  measurement  accuracy 
by  increasing  band  edge  data  while  decreasing  processing  time.  However,  the  direct 
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time-domain  system  may  need  additional  hardware  to  achieve  accurate  high  frequency 
results. 
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Figure  4.2:  Permittivity  of  MRAM  sample  using  S-band  rectangular  waveg¬ 
uide  for  the  forward  NRW  and  two  transmission  method. 
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Figure  4.3:  Permeability  of  MRAM  sample  using  X-band  rectangular  waveg¬ 
uide  for  the  forward  NRW  and  two  transmission  method. 


41 


Figure  4.4:  Permittivity  of  MRAM  sample  using  X-band  rectangular  waveg¬ 
uide  for  the  forward  NRW  and  two  transmission  method. 
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Figure  4.5:  Permeability  of  MRAM  sample  using  X-band  rectangular  waveg¬ 
uide  for  the  forward  NRW  and  two  transmission  method. 
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Figure  4.6:  Permittivity  of  MRAM  sample  measured  with  the  focus  arch 

free-space  system  from  2-6  GHz  for  the  forward  NRW  and  two  transmission 
method. 
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Figure  4.7:  Permeability  of  MRAM  sample  measured  with  the  focus  arch 
free-space  system  from  2-6  GHz  for  the  forward  NRW  and  two  transmission 
method. 
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Figure  4.8:  Permittivity  of  MRAM  sample  measured  with  the  focus  arch 

free-space  system  from  6-18  GHz  for  the  forward  NRW  and  two  transmission 
method. 
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Figure  4.9:  Permeability  of  MRAM  sample  measured  with  the  focus  arch 
free-space  system  from  6-18  GHz  for  the  forward  NRW  and  two  transmission 
method. 
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Figure  4.10:  Permittivity  of  MRAM  sample  measured  with  the  focus  arch 
free-space  system  from  2-6  GHz  for  the  frequency-domain  two  transmission 
method  and  time-domain  two  transmission  method. 
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Figure  4.11:  Permeability  of  MRAM  sample  measured  with  the  focus  arch 
free-space  system  from  2-6  GHz  for  the  frequency-domain  two  transmission 
method  and  time-domain  two  transmission  method. 
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Figure  4.12:  Permittivity  of  MRAM  sample  measured  with  the  focus  arch 
free-space  system  from  6-18  GHz  for  the  frequency-domain  two  transmission 
method  and  time-domain  two  transmission  method. 
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Figure  4.13:  Permeability  of  MRAM  sample  measured  with  the  focus  arch 
free-space  system  from  6-18  GHz  for  the  frequency-domain  two  transmission 
method  and  time-domain  two  transmission  method. 
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V.  Conclusions 


A  two  transmission  technique  was  developed  to  improve  constitutive  parameter  mea¬ 
surements.  The  permittivity  and  permeability  was  extracted  using  the  well  known 
NRW  algorithm  and  are  compared  to  the  two  transmission  technique  that  extracts 
permittivity  and  permeability  through  the  use  of  a  2D  Newton  root  search. 

The  technique  was  first  experimentally  demonstrated  using  a  MRAM  sample  in 
a  rectangular  waveguide  system  and  compared  to  the  NRW  method.  The  comparisons 
demonstrated  the  validity  of  this  method  and  the  improvement  over  the  NRW  method. 
The  technique  was  then  extended  to  a  frequency-domain  and  time-domain  focus  arch 
free-space  system.  A  differential  uncertainty  analysis  showed  that  due  to  the  path 
length  dependence  of  the  reflection  measurement  the  NRW  method  leads  to  greater 
uncertainty.  The  validity  of  using  a  time-domain  focus  arch  system  to  increase  band 
edge  data  was  demonstrated  and  compared  to  the  frequency-domain  system. 

The  results  of  this  thesis  show  the  improved  measurement  accuracy  of  consti¬ 
tutive  parameters  using  the  two  transmission  method  over  the  NRW  method.  Also, 
the  use  of  this  technique  in  a  time-domain  free-space  system  slightly  improves  band 
edge  data  while  removing  two  processing  operations. 

5.1  Future  Research 

There  are  a  few  general  improvements  that  could  be  made  to  increase  this 
techniques  application  to  direct  time-domain  measurements. 

5.1.1  Air  Gap  Error  Analysis.  A  further  investigation  into  the  error  due  to 
a  finite  air  gap  between  the  sample  and  acrylic  should  be  conducted.  This  analysis 
may  show  the  need  to  ensure  intimate  contact  between  the  sample  and  acrylic.  A 
three-layer  analysis  using  a  Newton  2D  root  search  can  be  conducted  similar  to  the 
two-layer  development  of  this  thesis. 

5.1.2  Low  Frequency  Application.  The  direct  time-domain  system  has  great 
potential  at  low  frequencies  in  its  current  configuration.  Low  frequency  preliminary 
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measurements  have  been  made  using  a  stripline  system.  Initial  results  show  improved 
band  edge  data  compared  to  the  conventional  NWA  frequency-domain  system.  The 
use  of  the  oscilloscope  and  TDT  module  should  produce  improved  low  frequency 
material  measurements. 

5.1.3  High  Frequency  Application.  The  direct  time-domain  system  may  be 
used  for  high  frequency  applications  with  additional  hardware.  The  use  of  a  pulse  gen¬ 
erator  and  remote  pulse  head  can  increase  the  rise  time  of  the  step  input.  Increasing 
the  rise  time  will  increase  spectral  content  and  enable  higher  frequency  applications. 
High  frequency  differential  time-domain  analysis  can  also  be  accomplished  using  Pi¬ 
cosecond  Pulse  Lab’s  TDR/TDT  module.  The  TDR/TDT  module  increases  the  40 
ps  rise  time  of  the  step  input  to  9  ps  resulting  in  much  greater  spectral  content. 
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Appendix  A.  ERAM  Results 


A.l  Waveguide  System 

A. 1.1  S-band  (2.6-3.95  GHz). 

A. 1.2  S-band  (8.2-12-4  GHz). 

A. 2  Waveguide  System 

A.  2.1  Frequency-Domain. 

A. 2. 2  Time-Domain. 
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Figure  A.l:  Complex  permittivity  of  125  mil  ERAM  sample  measured  with 
the  S-band  rectangular  waveguide  system  from  2.6-3.95  GHz  for  the  forward 
NRW  and  two  transmission  method. 
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Figure  A. 2:  Complex  permittivity  of  125  mil  ER.AM  sample  measured  with 
the  X-band  rectangular  waveguide  system  from  8.2-12.4  GHz  for  the  forward 
NRW  and  two  transmission  method. 
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Figure  A. 3:  Complex  permittivity  of  125  mil  ERAM  sample  measured  with 
the  focus  arch  free-space  system  from  2-6  GHz  for  the  forward  NRW  and  two 
transmission  method. 


57 


Figure  A. 4:  Complex  permittivity  of  125  mil  ERAM  sample  measured  with 
the  focus  arch  free-space  system  from  6-18  GHz  for  the  forward  NRW  and  two 
transmission  method. 
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Figure  A. 5:  Complex  permittivity  of  125  mil  ERAM  sample  measured 

with  the  focus  arch  free-space  system  from  2-6  GHz  for  the  frequency-domain 
method  and  direct  time-domain  method. 
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Figure  A. 6:  Complex  permittivity  of  125  mil  ERAM  sample  measured  with 
the  focus  arch  free-space  system  from  6-18  GHz  for  the  frequency-domain 
method  and  direct  time-domain  method. 
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